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Abstract 

We compute the one-loop QCD amplitudes for the decay of an off-shell vector boson 
with vector couplings into a quark-antiquark pair accompanied by two gluons keeping, 
for the first time, all orders in the number of colours. Together with previous work 
this completes the calculation of the necessary one-loop amplitudes needed for the 
calculation of the next-to-leading order 0(a^) corrections to four jet production in 
electron-positron annihilation, the production of a gauge boson accompanied by two 
jets in hadron-hadron collisions and three jet production in deep inelastic scattering. 





Multi-jet events in electron-positron annihilation have long been a source of vital infor¬ 
mation about the structure of QCD. In 1979, three jet events at DESY gave the first clear 
indications of the existence of the gluon [l|, while more recently, the colour factors that 
determine the gauge group of QCD have been determined [Q]. While the general features of 
events have been well described by leading order predictions, many quantities are known to 
suffer large radiative corrections. For example, the experimental four jet rate measured as a 
function of the jet resolution parameter y cut , was significantly higher than the leading order 
expectation. This is of particular interest at LEP II energies where 4 jet events are the main 
background to WW events where both W bosons decay hadronically. 


Recently there has been much progress towards a more complete theoretical description 
of four jet events produced in electron-positron annihilation. Dixon and Signer 0 have 
created a general purpose Monte Carlo program to evaluate the next-to-leading order cor¬ 
rections to a wide variety of four-jet and related observables. Such numerical programs have 
three important ingredients; the one-loop four parton amplitudes, the tree-level five parton 
amplitudes and a method for combining the four and five parton final states together. The 
tree-level amplitudes for e + e _ —> qqggg and e + e“ —> qqQQg are well known 0 , |,1 and sev¬ 
eral general methods have been developed for numerically isolating the infrared divergences 
0 !•!• However, the one-loop amplitudes are not fully known, and the next-to-leading 
order corrections described by the Dixon-Signer program 0] are therefore incomplete. 

The one-loop amplitudes relevant for this process are e + e~ —> qqgg and e + e“ —> qqQQ. 
Two groups BO have computed the amplitudes for the four quark final states keeping 
all orders in the number of colours, and these matrix elements have been implemented in 
the Dixon-Signer Monte Carlo [j3J. The situation for the two quark two gluon process is 
less complete and may be illustrated by examining the colour structure of the one-loop 
contribution to the four-jet cross section, 


<4gg° P = N 2 (N 2 - 1) (a< a) + 


—a® + —a® 
N 2 N 4 


( 1 ) 


where N is the number of colours. Bern, Dixon and Kosower II have provided compact 
analytic results for the leading colour part of the two quark two gluon process (<% a Q which 
should account for approximately 90% of the next-to-leading order corrections. Indeed, 
keeping all of the known contributions, Dixon and Signer find that the next-to-leading order 
predictions match onto the experimental data much better 0. However, the subleading 
terms a ® and ch c ) are still expected to contribute a further 10% to the 4-jet rate. In this 
letter, we compute the one-loop matrix elements for the 7 * —> qqgg process keeping, for the 
first time, all orders in the number of colours. These results can be implemented directly 
into the Dixon-Signer Monte Carlo [;fj, or can be used as a cross check with the results 
anticipated from the helicity approach |T3 . 


As in the calculation of the one-loop four quark matrix elements 13 and in contrast to 
the helicity method employed by Bern et al. [ 
i.e. the interference between tree-level and one-loop amplitudes. 


we compute ‘squared’ matrix elements, 
By doing so, we reduce all 
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tensor integrals down to scalars directly leaving at worst three (tensor) loop momenta in the 
n um erator of the box integrals. 


The evaluation of the tensor integrals represents the bulk of the calculation and should be 
treated with some care. The ‘squared’ tree-level matrix element for a given process contains 
inverse powers of certain invariants, resulting from the propagators present. This singularity 
structure is largely responsible for the infrared behaviour of the matrix elements. Based on 
general arguments, one expects that the one-loop corrections should have the same singu¬ 
larity structure as tree-level and should not contain more inverse powers of these (or other) 
invariants or additional kinematic singularities resulting from the calculational scheme used. 
However, as well as generating extra inverse powers of the invariants (above and beyond 
the expected tree-level singularities), traditional methods of reducing the tensor integrals to 
combinations of scalar integrals introduce (many) factors of the Gram determinant (A) in 
the denominator. Although the matrix elements themselves are finite as A —» 0, individual 
terms appear to diverge. These apparent singularities can be avoided by forming suitable 
combinations of the scalar integrals which are finite in the limit of vanishing A as well as in 
the limit where the invariant mass of a pair of particles tends to zero. Motivated by the work 
of Bern et al, [[L4|, we have shown in ref. |L5| how these functions are naturally obtained 
by differentiation of the original scalar loop-integral. The advantages of this decomposition 
are then two-fold. First, the expressions for the one-loop matrix elements are significantly 
simplified by grouping integrals together. Second, we can ensure that the one-loop singu¬ 
larity structure reproduces that at tree level. In other words, the only allowed kinematic 
singularities are those appearing at tree level, multiplied by coefficients that are well be¬ 
haved in all kinematic limits. This latter point is achieved by using identities amongst the 
set of basic functions and tends not to reduce the size of the answer, although the absence 
of further kinematic singularities ensures that the matrix elements are numerically stable. 
There is a slight price to pay in that, compared to the raw scalar integrals, the the set of 
basic functions is now enlarged and is not linearly independent, resulting in some ambiguity 
in the form of the final result. To generate and simplify our results, we have made repeated 
use of the algebraic manipulation packages FORM and MAPLE. 


For the process under consideration, 7 * —> qqgg we label the momenta as, 


7*Oi234) -► q(pi) + q(p2) + dips) + gfa), 


( 2 ) 


and systematically eliminate the photon momentum in favour of the four massless parton 
momenta. 

The colour structure of the matrix element at tree-level (n = 0) and one-loop (n = 1) is 
rather simple and we have, 


M (n) = e'WAM^L, = eg] ( —) 


- ILV3IP4 


2 n 


4:77/ 


x {(T“»T“) CICa 4 n) (3,4) + (T‘“T‘=) cta A^(i, 3) + ±6 cia 6 am A ( p(3, 4) J , 


(3) 
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where Ci, C 2 are the colours of the quarks and 03 , CI 4 the colours of the gluons. The arguments 
of Ai indicate a permutation of the momenta of the external gluons. At lowest order, 

4 0 ) (3,4) = 0, (4) 

while at one-loop we find, 

■Af’M = iV.4 1 , (3,4)-i4 1 ) (3,4), (5) 

= >(g>(3,4)=>(g ) (4,3). ( 6 ) 

The functions a = A, B,C represent the contributions of the three gauge invariant 

sets of Feynman diagrams shown in Figs. 1, 2 and 3 respectively. We note that the contri¬ 
bution 4i (3,4) also includes terms proportional to rip/N, where np is the number of light 
quark flavours. These terms originate from the fermionic self-energy and triple gluon vertex 
corrections shown in the last two diagrams of Fig. 1. The Dixon-Signer Monte Carlo [}3[] 
incorporates the leading colour term, M^(3,4). 


The squared matrix elements are relatively straightforwardly obtained at leading order 
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Some care must be taken to ensure that only physical, transverse degrees of freedom 
are included for the external gluons. We do this by requiring that terms in that 

vanish when contracted with the physical gluon polarisation vectors e 1 ' 3 or e U4 are removed 
by hand [|T7| . The sum over the polarisations of all the gauge bosons may then be performed 
in the Feynman gauge, 


E 

spins 


e'V 1 ' 


= -<r ■ 


(7) 


ffence, 

E l^ (0) I 2 = C ^(N 2 - l)iv(T(3, 4) + T(4, 3) - E T | 

spins l ) 

where, 

T(3,4)= £ |4 0)t (3,4)4 0) (3,4)|, 

spins 

and, 

T= £ |(Af )t (3,4) + 4° )t (4,3)) (4 0) (3,4) + 4 0) (4,3)) 

spins 


( 8 ) 

(9) 

( 10 ) 


The 3-gluon vertex contributions to (3,4) and 4°^(4, 3) enter with opposite sign, so T 
(with no arguments) is the contribution from the pure QED-likc diagrams. 


The relevant ‘squared’ matrix elements are the interference between the tree-level and 
one-loop amplitudes, 

5] ^ ( N 2 - 1 )N 

spins 

x j£ A (3,4) + C A { 4,3) - E (c A + C B (3, 4) + C B ( 4,3) - C c ) + ^ fl |, (11) 
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(b) 


Figure 1: The classes of Feynman diagrams relevant for the function *4.^(3,4). Reading 
clockwise round the diagram and starting from the quark (pi) at the top, we encounter 
gluon (p 3 ) before gluon (p 4 ) and end at antiquark ( p 2 ). The solid circle indicates the possible 
positions for attaching the off-shell photon to the quark-antiquark pair. Diagrams (a), taken 
with both permutations of gluons 3 and 4, contribute to the piece Ca while the permutation 
shown in (a) + (b) gives the contribution to £a( 3,4). Diagrams with self-energy corrections 
on the external lines are zero in dimensional regularisation and have been omitted. 
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Figure 2: The classes of Feynman diagrams relevant for the function *4.^(3,4) Reading 
clockwise round the diagram and starting from the quark (pi) at the top, we encounter 
gluon (p 3 ) before gluon (p 4 ) and end at antiquark (p 2 ). The solid circle indicates the possible 
positions for attaching the off-shell photon to the quark-antiquark pair. Diagrams (a), taken 
with both permutations of gluons 3 and 4, contribute to the piece C B while the permutation 
shown in (a) + (b) gives the contribution to £#(3,4). 
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Figure 3: The classes of Feynman diagrams relevant for the function Wq when taken with 
both permutations of the gluons. Reading clockwise round the diagram and starting from 
the quark {p\) at the top, we encounter gluon (p 3 ) before gluon (p 4 ) and end at antiquark 
(p 2 )- The solid circle indicates the possible positions for attaching the off-shell photon to 
the quark-antiquark pair. 


with, 

£«(3,4)= £ |4 1)t (3,4)4 0) (3,4)|, 

spins 

for a = A, B and the QED-like structures, 

£«= E IK" (3.4)+^* (4,3)) (<’(3,4) + <'(4,3)) 

spins 


( 12 ) 


(13) 


Hence the ‘squared’ matrix elements are described by 5 independent jC a . The £ a (3,4) are 
symmetric under the exchange p\ p 2 and p 3 p A while the C a are symmetric under 
either p\ p 2 or p 3 p A . 


Working in dimensional regularisation with 4 — 2e spacetime dimensions, it is straight¬ 
forward to remove the infrared and ultraviolet poles from the C a since they are proportional 
to the tree-level amplitudes, 


£.4(3,4) 

C A 


f £3 £*24 'P 34 

V £ 7 £ 

1 / £*13 ^24 

2 £ _ 

1 / V\4 V 23 

2 V £ £ ~ 


£3,4) + £.4(3,4), 

5 1 _ l 1 - 2 ? 2 ) (T+T(3,4) _ T(4,3)) 

~ £ + £ 4 > 3 ) - £3,4)) + C A , 


(14) 

(15) 
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C B { 3,4) = 
£b = 
£c = 


V 


12 


T(3,4)+£b(3,4), 


P 


12 


P 


12 


P 1234 ^ 

)T+£ b , 


e v 

) 

P 34 

, V 14 , 

P 23 

" I 2 ” 

+ — + 
e z 

l 2- 

P 34 

, P 13 . 

P 24 


(16) 

(17) 


(T + T(4,3)-T(3,4))+£ c . (18) 


where we have introduced the notation, 

'47r/i 2 \ £ r 2 (l-e)r(l + e) 


V,= 


1(1 - 26) 


(19) 


In physical cross sections, this pole structure cancels with the infrared poles from the 7 * —> 
qqgg + g process and those generated by ultraviolet renormalisation. 

In determining the finite pieces, C, we are concerned to ensure that the singularity 
structure matches that of the tree-level functions T(3,4) and T. I 11 terms of the generalised 
Mandelstam invariants, 


Sij = G Vi + Pj), Sij k = (Pi + pj + Pk) , s ijki = G Pi + Pj +Pk + Pi) , 


( 20 ) 


T contains single poles in S13, S23, S14 and S24 while T(3,4) has poles in S13, S34 and S24. In 
addition, both T and T (3,4) contain double poles in the triple invariants S 134 and S 234 . Using 
the tensor reduction described in |0J, possible singularities due to Gram determinants are 
automatically protected. However, it is possible to generate apparent singularities in double 
or triple invariants such as Si 2 or S 123 . These poles do not correspond to any of the allowed 
infrared singularities and the matrix elements are finite as, for example, S 12 —> 0 or S 123 —*► 0 . 
In fact, it is straightforward to explicitly remove such poles using identities amongst the 
combinations of scalar integrals. For example, the identity, 


(^123 — ^12) 


Sl23 


Ld 22 (pi2,p3,p4) = -Ld 2 l(pi2,P3,P4) + Ldi(pi2,p 3 ,p4) + Lei (pi2, p 3 4.) , (21) 


relating the functions for box integral functions with two adjacent masses (defined in Ap¬ 
pendix B of 0) is useful to eliminate poles in S 123 . The finite pieces can be written 
symbolically as, 

£ = E^(s)L<, (22) 

i 

where the coefficients T)(s) are rational polynomials of invariants. The finite functions 
L i are the linear combinations of scalar integrals defined in JT5| which are well-behaved 
in all kinematic limits. Any denominators of the corresponding tree-level matrix element 
are allowed in the coefficient Pj(s), with any additional fake singularities protected by L*. 
Typically the coefficient of a given function contains 0(100) terms, comparable with the size 
of the tree-level matrix elements. The number of functions is rather large, of 0(100), but 
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is a minimal set which protects all the kinematic limits and is therefore numerically stable. 
The analytic expressions for the individual C a , and particularly the subleading colour terms 
Cb and Cc, are therefore rather lengthy so we have constructed a FORTRAN subroutine to 
evaluate them for a given phase space point. This can be directly implemented in general 
purpose next-to-leading order Monte Carlo programs for the e + e - —> 4 jet, e ± p —> + 3 jet 
and pp —> V + 2 jet processes. 


To summarize, we have performed the first calculation of the one-loop ‘squared’ matrix 
elements for the 7 * —> qqgg process keeping all orders in the number of colours. We have 
grouped the Feynman diagrams according to the colour structure, but unlike the helicity ap¬ 
proach of [|T|], have used conventional dimensional regularisation throughout. We have used 
the one-loop reduction method described in ||T5| to obtain the finite parts of the one-loop 
matrix amplitudes in terms of functions that are well behaved in all kinematic limits. The 
finite one-loop expressions C a contain the same singularity structure as the tree level ampli¬ 
tudes T and are numerically stable. Our results are rather lengthy and we have provided a 
FORTRAN implementation of them that can be either used with the existing Dixon-Signer 
program for e + e _ —> 4 jets or as a completely independent check of amplitudes obtained 
using the helicity approach |[I], [12], [L3|. Together with previous work BB0 this com¬ 
pletes the calculation of the necessary one-loop amplitudes for the coupling of an electroweak 
gauge boson to four massless partons with the following exceptions: 


(a) Contributions proportional to the axial coupling. 


(b) Contributions where the electroweak boson couples to a closed fermion loop. 

Both of these contributions are expected to be small because of cancellations between the 
up- and down-type quark contributions. 
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